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ABSTRACT 



We present NLTE Li abundances for 88 stars in the metallicity range -3.5 < [Fe/H] < -1.0. The effective temperatures are based on the 
infrared flux method with improved E(B - V) values obtained mostly from interstellar NaiD lines. The Li abundances were derived 
through MARCS models and high-quality UVES+VLT, HIRES-l-Keck and FIES+NOT spectra, and complemented with reliable 
equivalent widths from the literature. The less-depleted stars with [Fe/H]< -2.5 and [Fe/H] > -2.5 fall into two well-defined plateaus 
of Au = 2.18 (cr = 0.04) and Au = 2.27 (a = 0.05), respectively. We show that the two plateaus are flat, unlike previous claims for 
a steep monotonic decrease in Li abundances with decreasing metallicities. At all metallicities we uncover a fine-structure in the Li 
abundances of Spite plateau stars, which we trace to Li depletion that depends on both metallicity and mass. Models including atomic 
diffusion and turbulent mixing seem to reproduce the observed Li depletion assuming a primordial Li abundance Au = 2.64, which 
agrees well with current predictions (Au = 2.72) from standard Big Bang nucleosynthesis. Adopting the Kurucz overshooting model 
atmospheres increases the Li abundance by H-0.08 dex to Au = 2.72, which perfectly agrees with BBN-l-WMAP. 
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1. Introduction 

One of the most important discoveries in the study of the chem- 
ical composition of stars was made in 1982 by Monique and 
Francois Spite, who found an essentially constant Li abundance 
^-_l in warm metal-poor stars (Spite & Spite 1982), a result inter- 

^ preted as a relic of primordial nucleosynthesis. Due to its cos- 
■"sj" mological significance, there have been many studies devoted to 
■^ Li in metal-poor field stars (e.g., Ryan et al. 1999; Melendez & 
0\ Ramirez 2004, hereafter MR04; Boesgaard et al. 2005, hereafter 
^. BOS; Charbonnel & Primas 2005; Asplund et al. 2006, hereafter 
l/-^ A06; Shi et al. 2007, hereafter S07; Bonifacio et al. 2007, here- 
O after B07; Hosford et al. 2009; Aoki et al. 2009; Sbordone et al. 
^D . 2010, hereafter SI 0), with observed Li abundances at the lowest 
^H . metallicities ([Fe/H] ~ -3) from Au = 1.94 (B07) to Au = 2.37 

'^ ; (MR04). 

^ A primordial Li abundance of Au = 2.72;^[![!^ is predicted 



X 



0.06 

(Cyburt et al. 2008; see also Steigman 2009; Coc & Vangioni 
2010) with the theory of big bang nucleosynthesis (BBN) and 
the baryon density obtained from WMAP data (Dunkley et al. 
2009), which is a factor of 2-6 times higher than the Li abun- 
dance inferred from halo stars. There have been many theoretical 
studies on non-standard BBN trying to explain the cosmological 
Li discrepancy by exploring the frontiers of new physics (e.g. 
Coc et al. 2009; locco et al. 2009; Jedamzik & Pospelov 2009; 
Kohri & Santoso 2009). Alternatively, the Li problem could be 
explained by a reduction of the original Li stellar abundance due 
to internal processes (i.e., by stellar depletion). In particular, stel- 



* Based in part on observations obtained at the W. M. Keck 
Observatory, the Nordic Optical Telescope on La Palma, and on data 
from the HIRES/Keck archive and the European Southern Observatory 
ESO/ST-ECF Science Archive Facility 



lar models including atomic diffusion and mixing can deplete a 
significant fraction of the initial Li content (Richard et al. 2005; 
Piau 2008), although such models depend on largely uncon- 
strained free parameters. On the other hand, it is not easy to rec- 
oncile the lack of observed abundance scatter in the Spite plateau 
with substantial Li depletion (e.g. Ryan et al. 1999; A06). Due to 
the uncertainties in the Li abundances and to the small samples 
available, only limited comparisons of models of Li depletion 
with stars in a broad range of mass and metallicities have been 
performed (e.g., Pinsonneault et al. 2002; B05). 

The observed scatter in derived Li abundances in previous 
studies of metal-poor stars can be as low as 0.03 dex (e.g. Ryan et 
al. 1999; A06), fully consistent with the expected observational 
errors. Yet, for faint stars observational errors as high as 0.2 dex 
have been reported (e.g., Aoki et al. 2009). In order to provide 
meaningful comparisons with stellar depletion models, precise 
Li abundances for a large sample of stars are needed. Here we 
present such a study for the first time for a large sample of metal- 
poor stars (-3.5 < [Fe/H] < -1.0) with masses in a relatively broad 
mass range (0.6-0.9 Mq). 



2. Stellar sample and atmospheric parameters 

The sample was initially based on the 62 stars analysed by 
MR04. We added stars with published equivalent width (EW) 
measurements based on high-quality observations (Sect. 3), and 
included also stars for which we could obtain EW Li measure- 
ments from UVESh-VLT, HIRES +Keck and FIES-)-NOT spectra 
(either from our own observations or from the archives). The re- 
solving power of the different instruments used ranges from ~ 
45,000 to 110,000, and the S/N is typically > 100. Some mea- 
surements taken from the literature were obtained with lower re- 
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solving power (-35,000), or high resolving power (R ~ 50,000) 
but lower S/N. When necessary we averaged several measure- 
ments in order to decrease the errors. The stars for which we 
could not obtain EW better than ~5% were discarded, with a few 
exceptions (see Sect. 3). Our adopted log g values and metallici- 
ties [Fe/H] are mean values (after discarding outliers) taken from 
the compilation of stellar parameters by Melendez (in prepara- 
tion), who has updated the metallicity catalogue of Cayrel de 
Strobel et al. (2001) from 4918 entries to more than 14000. We 
only keep stars with at least two previous spectroscopic anal- 
yses, so that their metallicities and surface gravities have been 
confirmed at least by one other study. The sample has 88 stars 
with stellar parameters 5250 < Teff < 6600 K, 3.6 <logg < 4.8 
and -3.5 < [Fe/H] < -1.0. 

We obtained effective temperatures (Teff) by a new im- 
plementation (Casagrande et al. 2010, hereafter CIO; see also 
Casagrande et al. 2006) of the infrared flux method (IRFM). In 
our previous work on Li abundances in metal-poor stars (MR04), 
Teff was obtained from the IRFM implementation by Ramirez & 
Melendez (2005), relying on the calibration of the bolometric 
correction by Alonso et al. (1995), which contains only a few 
low-metallicity stars; thus, for very metal-poor turn-off stars, it 
is not well defined and may introduce systematic errors. Indeed, 
CIO have found that below [Fe/H] ^ -2.5, the Ramirez & 
Melendez (2005) T^ff go from being too cool by ^ lOOK at 
[Fe/H]=-2.5 to too hot by ^ lOOK at [Fe/H]=-3.5. 

Our new IRFM implementation determines the bolometric 
flux from observed broad-band photometry (BVRIJHK), thus 
representing an important improvement over the T^ff obtained 
by Ramirez & Melendez (2005). Furthermore, the zero point of 
our Teff scale has been carefully checked for accuracy using so- 
lar twins (their average Teff matching the solar one within a few 
K), spectral energy distributions measured with HST for a metal- 
rich star with planet and a metal-poor turn-off star, and stars with 
interferometric angular diameters (for which our IRFM T^ff and 
fluxes result in angular diameters only 0.6 + 1.7% smaller, corre- 
sponding to 18+50 K). Thus, contrary to most previous studies, 
the uncertainty of our absolute Li abundances is not dominated 
by the unknown zero point of the adopted Tetf scale. 

In order to obtain Teff values with IRFM, one needs good 
estimates of E{B - V) to properly correct the photometric data. 
In our previous work we used extinction maps (see Sect. 4 in 
Melendez et al. 2006, hereafter M06), but in order to achieve 
a very high precision (~ 0.005 mag), we must rely on other 
methods. Thanks to the large radial velocity of most metal- 
poor stars, the interstellar NaiD lines can be used to obtain 
this precision. We employed this technique (Alves-Brito et al. 
2010) when high-resolution spectra were available (~72% of the 
cases), otherwise we used extinction maps (M06). Our adopted 
E(B - V) values are shown in Fig. 1 and have median errors of 
0.003+0.002 mags, implying errors in Teff of ±15 K. Adding in 
errors in photometry, bolometric flux, log g and [Fe/H], we esti- 
mate a typical total internal error in T^fi of +50 K. Q 

We obtained stellar masses from the a-enhanced Y^ 
isochrones (Demarque et al. 2004). We interpolated a fine grid 
of models with a step A[Fe/H] = 0.02 dex and adopting [a/Fe] 
= for [Fe/H] > 0, [a/Fe] = -0.3 x[Fe/H] for -1 < [Fe/H] < 
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Fig. 1. Adopted E(B-V) values vs. [Fe/H]. Filled squares are val- 
ues based on NaD lines and open squares on reddening maps. 



0, and [c/Fe] = +0.3 for [Fe/H] < -l.Q At a given metallicity, 
we searched for all solutions allowed by the error bars in T^ff, 
Hipparcos parallaxes (reliably available for 2/3 of the sample) 
and [Fe/H], adopting the median value in age and mass, and the 
standard deviation of the solutions is adopted as the error. We 
also used the mean log g values found in the literature to esti- 
mate another set of masses and ages. The masses obtained by 
both methods agree very well (median difference [Hipparcos - 
literature] of -0.001 Mq). The weighted averages of both results 
were adopted. We find typical ages of ~ 1 1 Gyr and masses in 
the range of 0.6-0.9 Mq. In a few cases the adopted stellar pa- 
rameters resulted in low ages (< 6 Gyr), probably due to errors 
in log g. In those cases we estimated masses adopting an age=l 1 
Gyr, but the choice of age does not affect our conclusions (e.g., 
Aage=+1 Gyr results in Amass=-0.007 M© for a dwarf with 
[Fe/H]=-2 and Teff «i6300 K). 



3. Li abundance analysis 

The LTE Li abundances were obtained with our code (A06) with 
MARCS (Gustafsson et al. 2008) models and the 2002 version 
of MOOG (Sneden 1973) with Kurucz overshooting (Castelli et 
al. 1997) models (which have improved overshooting with re- 
spect to earlier models that showed several problems described 
in van't Veer-Menneret & Megessier 1996). Both model atmo- 
spheres give similar Li abundances except for a zero-point differ- 
ence of ~0.08 dex (Kurucz overshooting - MARCS). According 
to MR04 the difference in Au between Kurucz overshooting 
and no-overshooting models is also ~0.08 dex, so Kurucz no- 
overshooting models give similar Au as the MARCS models. 
The results presented herein will be based on the MARCS mod- 
els, because the overshoot option implemented by Kurucz is not 
supported by 3D hydrodynamical model atmospheres (Asplund 
et al. 1999; see also Heiter et al. 2002); we will nevertheless 
mention the implication of using Kurucz overshooting models 
for completeness. We adopt the ID non-LTE corrections by Lind 
et al. (2009a), which are based on the most recent radiative and 
collisional data (Barklem et al. 2003). The difference between 
LTE and non-LTE is relatively constant for our stars so that sim- 
ilar conclusions would have been obtained assuming LTE, yet, 
for an accurate analysis NLTE must be taken into account be- 
cause it affects the absolute level of the Spite plateau. For metal- 
poor ([Fe/H] < -1) stars with Teff > 6000 K the typical NLTE 
correction is -0.06 + 0.01 dex. 



' In some cases the Teff given in Table 1 may be slightly different 
from the T^ff given in C 10 due to improved input parameters (e.g., better 
E(B-V) values) in the IRFM used here. The conclusions of CIO are 
unaffected by these minor changes. 



- The adopted relation between [ff/Fe] and [Fe/H] ignores the de- 
tailed abundance patterns of the two discs and halos (e.g., Ramirez et 
al. 2007, Nissen & Schuster 2010). 
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Fig. 2. Li abundances vs. reft for our sample stars in different 
metallicity ranges. A typical error bar is shown. 



We measured the EW of the Li feature at 6707.8 A us- 
ing UVES+VLT, HIRES+Keck and FIES+NOT high-resolution 
spectra, from which we obtain typical errors in EW of 0.7 mA, 
which includes uncertainties in continuum placement. We com- 
plemented our measurements with data from the literature (BOS; 
A06; B07; S07; Asplund & Melendez 2008, hereafter AM08; 
SIO). An important improvement with respect to our previous 
work is that now we select stars with errors in EW below 5% 
(typically ~ 2-3%), instead of the 10% Umit adopted in MR04. 
The only exceptions are the cool dwarfs HD64090 and BD-i-38 
4955, which are severely depleted in Li and have EW errors of 
8% and 10%, respectively, and the very metal-poor stars (B07) 
BPS CS29518-0020 (5.2%) and BPS CS29518-0043 (6.4%), 
which were kept due to their low metallicity. 

The main sources of error are the uncertainties in equivalent 
widths and T^ff, which in our work have typical values of only 
2.3% and 50 K, implying abundance errors of 0.010 dex and 
0.034 dex, respectively, and a total error in Au of ~ 0.035 dex. 
This low error in Au is confirmed by the star-to-star scatter of the 
Li plateau stars, which have similar low values (e.g. cr - 0.036 
dex for [Fe/H] < -2.5, see below). Our Li abundances and stellar 
parameters are given in Table 1 (online material). 



4. Discussion 

4.1. The Tf,s cutoff of the Spite plateau 

Despite the fact that Li depletion depends on mass (e.g. 
Pinsonneault et al. 1992), this variable has been ignored by most 
previous studies. Usually a cutoff in Teff is imposed to exclude 
severely Li-depleted stars in the Spite plateau, with a wide range 
of adopted cutoffs, such as 5500 K (Spite & Spite 1982), 5700 K 
(B05), 6000 K (MR04; S07) and -6100 K for stars with [Fe/H] 
< -2.5 (Hosford et al. 2009). 

At a given mass, the Teff of metal-poor stars has a strong 
metallicity-dependence (e.g. Demarque et al. 2004). As shown 
in Figs. 1 1-12 of M06, the Tgft of turnoff stars increases for de- 
creasing metallicities. Hence, a metallicity-independent cutoff in 
Jeff may be an inadequate way to exclude low-mass Li-depleted 
stars from the Spite plateau. As show in Fig. 2, where Au in dif- 
ferent metallicity bins is shown as a function of Tgft, stars with 
lower Teff in a given metallicity regime are typically the stars 
with the lowest Li abundances, an effect that can be seen even 

in the sample stars with the lowest metallicities ([Fe/H] 3). 

This is ultimately so because the coolest stars are typically the 
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Fig. 3. Li abundances for stars with T^g > 5700 K (open circles), 
> 6100 K (filled squai-es), > 6350 K (filled tiiangles) and > 5850 
- 180x[Fe/H] (stars). In the bottom panel stars above the cut- 
off in Teft fall into two flat plateaus with cr=0.04 and 0.05 dex 
for [Fe/H] < -2.5 (dotted line) and [Fe/H] > -2.5 (solid line), 
respectively. 



least massive, and therefore have been more depleted in Li (see 
Sect. 4.3). 

In Fig. 3 we show the Li abundance for cutoffs = 5700 K 
(open circles), 6100 K (filled squares) and 6350 K (filled trian- 
gles). Using a hotter cutoff is useful to eliminate the most Li- 
depleted stars at low metallicities, but it removes from the Spite 
plateau stars with [Fe/H] > -2. Imposing a hotter Teft cutoff at 
low metallicities and a cooler cutoff at high metallicities elim- 
inates the most Li-depleted stars at low metallicities, but keeps 
the most metal-rich stars in the Spite plateau. We propose such 
a metallicity-dependent cutoff below. 

4.2. Two flat Spite plateaus 

Giving the shortcomings of a constant Tetf cutoff, we propose an 
empirical cutoff of Tgft = 5850 - 180x[Fe/H]. The stars above 
this cutoff are shown as stars in the bottom panel of Fig. 3. 
Our empirical cutoff excludes only the most severely Li-depleted 
stars, i.e., the stars that remain in the Spite plateau may still be 
affected by depletion. The less Li-depleted stars in the bottom 
panel of Fig. 3 show two well-defined groups separated at [Fe/H] 
2.5 (as shown below, this break represents a real discontinu- 
ity), which have essentially zero slopes (within the error bars) 
and very low star-to-star scatter in their Li abundances. The first 
group has -2.5 < [Fe/H] < -1.0 and <ALi>i = 2.272 (o-=0.051) 
dex and a slope of 0.018+0.026, i.e., flat within the uncertainties. 
The second group is more metal poor ([Fe/H] < -2.5) and has 
<ALi>2 - 2.184 dex (cr = 0.036) dex. The slope of this second 
group is also zero (-0.008±0.037). Adopting a more conserva- 
tive exponential cutoff obtained from Y^ isochrones (Demarque 
et al. 2004), which for a 0.79 Mq star can be fit by T^s = 6698 
-2173 X elF<=/H]/i.02i^ .^g would also recover a flat Spite plateau. 
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Fig. 4. (a) filled circles show the slope/cr of Ali vs. [Fe/H] in the 
range Xmin < [Fe/H] < -1.0, where Xmin is in the interval [-3.5,- 
1.25]. For all values of X^in > -2.5, the slope is insignificant 
(« 3cr); only when stars with [Fe/H] < -2.5 are included, a 
slope is forced in Au vs. [Fe/H]. Open circles show the slope/o" 
for the range -3.25 < [Fe/H] < X^ax, where -3.0 < X^ax ^ 
-1.0. ForXmax < -2.5 the slope is negligible (<< 3cr), but when 
stars with [Fe/H] > -2.5 are included, a slope is generated, (b) 
The correlation coefficient rspearman is shown by triangles, with 
filled and open symbols with similar meaning as in panel (a). The 
probabilities of a correlation between Ali and [Fe/H] by pure 
chance, associated to the filled and open triangles, are shown by 
solid and dashed lines, respectively. No meaningful correlation 
(probability « 1) exists for either [Fe/H] < -2.5 or [Fe/H] > 
-2.5. 



although only stars with [Fe/H] > -2.5 are left using this more 
restrictive cut-off. Thus, the flatness of the Spite plateau is inde- 
pendent of applying a linear or an exponential cutoff. 

Adopting a constant cutoff in Teff we also find flat plateaus. 
For example adopting a cutoff of T^s > 6100 K (filled squares in 
Fig. 3) we find in the most metal-rich plateau ([Fe/H] > -2.5) no 
trend between Li and [Fe/H] (slope - 0.019±0.025, Spearman 
rank correlation coefficient rspearman=0.1 and a probability of 
0.48 (i.e., 48%) of a correlation arising by pure chance for 
[Fe/H] > -2.5), while for the most metal-poor plateau ([Fe/H] 
< -2.5) we also do not find any trend within the errors (slope 
= 0.058+0.072, rspeaiman=0.2 and 41% probability of a spurious 
correlation). Using a hotter cutoff {T^s > 6350 K, filled trian- 
gles in Fig. 3) we obtain also two flat plateaus with slope - 
-0.040+0.063 (rspearman = -0.2, probability = 60%) for [Fe/H] 

> -2.5 and slope = 0.008+0.035 (rspearman = 0.1, probability = 
68%) for [Fe/H] < -2.5. 

In Fig. 4 we demonstrate that the break at [Fe/H] ~ -2.5 
is statistically significant. In panel (a) the slope/cr in the Ali vs. 
[Fe/H] plot for the range X^in < [Fe/H] < -1.0 are shown as 
filled circles, where Xmi„ varies within [-3.50,-1.25]. For X^in 

> -2.5, the slope is insignificant (<< 3cr), and only when stars 
with [Fe/H] < -2.5 are included a measurable slope is forced in 
the Au vs. [Fe/H] relation. The opposite test is shown by open 
circles, where we show slope/cr for the range -3.25 < [Fe/H] < 



Xmax, where Xtnax changes from [-3.00,-1.00].ForXn,ax < -2.5 
the slope is negligible (<< 3cr), and only when stars with [Fe/H] 
> -2.5 are included, a slope is produced between Au and [Fe/H]. 
The correlation coefficient rspearman and the probability of a cor- 
relation between Ali and [Fe/H] by pure chance are shown in 
panel (b). Again, this plot shows that no correlation between Ali 
and [Fe/H] exists within the two groups (-3.5 < [Fe/H] < -2.5 
and -2.5 > [Fe/H] > -1.0), and that only when stars between the 
two groups are mixed, significant (probability ~ 0) correlations 
of Au with [Fe/H] are generated. Systematically lower E(B-V) 
values (by ~0.03 mags) in the most metal-poor plateau ([Fe/H] 
< -2.5) could produce a more Li-depleted plateau. However, 
those E(B-V) values (Fig. 1) are actually slightly higher (~0.007 
mags, i.e., ~33 K) than those for [Fe/H] > -2.5 (mainly due to a 
high number of unreddened nearby more metal-rich stars), thus 
not explaining the existence of two different plateaus. 

Previous claims of a steep monotonic decrease in Li abun- 
dance with decreasing metallicity (e.g. Ryan et al. 1999; A06; 
B07; Hosford et al. 2009) are probably due to the mix of stars 
from the two different groups, forcing a monotonic dependence 
with metallicity. Our large sample of homogeneous and precise 
Li abundances that covers a broad metallicity range (-3.5 < 
[Fe/H] < -1.0) does not support these claims. Nevertheless, a 
hint of two different groups in the Spite plateau was already 
found by A06, who found a change in the slope of the Spite 
plateau at [Fe/H] ^ -2.2. Also, in the combined A06+B07 sam- 
ple (Fig. 7 of B07), there are two different groups: stars with 
[Fe/H] > -2.6 have Au > 2.2, while stars with [Fe/H] < -2.6 
have Au < 2.2. Although in the study by MR04 a flat Spite 
plateau is found in the range -3.4 < [Fe/H] < -1, this is due 
to the overestimation of Teff below [Fe/H] < -2.5, thus overes- 
timating Ali at low metallicities and forcing a flat plateau from 
[Fe/H] = -3.4 to -1. 



4.3. Correlation between Li and mass 

Models of Li depletion predict that the least massive stars are 
the most depleted in Li, but due to the limitations of previous 
samples, these predictions have not been thoroughly tested at 
different metallicity regimes in metal-poor starsQ In Fig. 5 we 
show our Li abundances as a function of stellar mass for dif- 
ferent metallicity ranges. As can be seen, the Li plateau stars 
have a clear dependence with mass for all metallicity regimes. 
Excluding stars with mass < O.7M0 (including those stars will 
result in even stronger correlations), linear fits result in slopes 
of 6, 3, 2, 2 dex Mg' for stars in the metallicity ranges -1.3 < 
[Fe/H] < -1.0, -1.6 < [Fe/H] < -1.3, -2.5 < [Fe/H] < -1.6, 
and -3.5 < [Fe/H] < -2.5. The slopes are significant at the 8, 2, 
5, 1 cr level, respectively. The correlation coefficient rspearman is 
0.9, 0.6, 0.6, 0.3, and the probability of a correlation between Ali 
and mass arising by pure chance is very small: 5x10"^, 3x10"^, 
IxlO""*, and 1.3x10"', for stars in the same metallicity ranges as 
above. Thus, the correlations of Ali and mass in different metal- 
licity regimes are very significant. 

Recently, Gonzalez Hernandez et al. (2008, hereafter G08) 
have studied the metal-poor ([Fe/H] 3.5) double-lined spec- 
troscopic binary BPS CS22876-032, providing thus crucial data 
to test our Li-mass trend. With our method and the stellar param- 
eters of G08, we obtain a mass ratio of 0.89, very close to their 
value (0.911+0.022) obtained from an orbital solution. For the 



^ Except for the work of BOS, who provide comparisons at different 
metallicities but for cool severely Li-depleted dwarfs, i.e. probably of 
lower mass than most stars shown in Fig. 5 
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primary we obtain Ma - 0.776 Mq, and adopting the mass ratio 
of G08, Mb = 0.707 Mq is obtained for the secondary. The LTE 
Li abundances were taken from G08 and corrected for NLTE ef- 
fects ( 0.05 dex). The components of the binary are shown as 

circles in Fig. 5, nicely following the trend of the most metal- 
poor stars. Including this binary in our sample would strengthen 
the Li-mass correlation of stars with [Fe/H] < -2.5. A slope = 3 
dex Mq' significant at the 3cr level is obtained, with rspearman - 
0.5 and a low probability (2x10^^) of the trend being spurious. 

While mass-dependent Li depletion is expected from stan- 
dard models of stellar evolution, this should only occur at signif- 
icantly lower masses than considered here. These stellar models 
only predict very minor ^Li depletion (^ 0.02 dex) for metal- 
poor turn-off stars (e.g. Pinsonneault et al. 1992), which is far 
from sufficient to explain the ~0.5dex discrepancy between 
the observed Li abundance and predictions from BBN+WMAP. 
Bridging this gap would thus require invoking stellar models that 
include additional processes normally not accounted for, such as 
rotationally-induced mixing or diffusion. In Fig. 5 we confront 
the predictions of Richard et al. (2005) with our inferred stel- 
lar masses and Li abundances. The models include the effects 
of atomic diffusion, radiative acceleration, and gravitational set- 
tling, but moderated by a parameterized turbulent mixing (T6.0, 
T6.09, and T6.25, where higher numbers mean higher turbu- 
lence); so far only the predictions for [Fe/H] =-2.3 are available 
for different turbulent mixing models. The agreement is very 
good when adopting a turbulent model of T6.25 (see Richard 
et al. for the meaning of this notation) and an initial Au - 2.64; 
had the Kurucz convective overshooting models been adopted, 
the required initial abundance to explain our observational data 
would correspond to Au = 2.72. Two other weaker turbulence 
models that produce smaller overall Li depletions are also shown 
in Fig. 5, but they are less successful in reproducing the observed 
Li abundance pattern. Our best-fitting turbulence model (T6.25) 
is different from that (T6.0) required to explain the abundance 
pattern in the globular cluster NGC 6397 at similar metallicity 
(Korn et al. 2006; Lind et al. 2009b). Another problem with 
adopting this high turbulence is that the expected correspond- 
ing ^Li depletion would amount to > 1.6 dex and thus imply 
an initial ^Li abundance at least as high as the primordial ^Li 
abundance if the claimed ^Li detections in some halo stars (A06; 
AMOS) are real. 

Our results imply that the Li abundances observed in Li 
plateau stars have been depleted from their original values and 
therefore do not represent the primordial Li abundance. It ap- 
pears that the observed Li abundances in metal-poor stars can be 
reasonably well reconciled with the predictions from standard 
Big Bang nucleosynthesis (e.g. Cyburt et al. 2008) by means of 
stellar evolution models that include Li depletion through dif- 
fusion and turbulent mixing (Richard et al. 2005). We caution 
however that although encouraging, our results should not be 
viewed as proof of the Richard et al. models until the free pa- 
rameters required for the stellar modelling are better understood 
from physical principles. 
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Fig. 5. Li abundances as a function of stellar mass in different 

metallicity ranges. The metal-poor ([Fe/H] 3.5) binary BPS 

CS22876-032 (GOB) is represented by open circles, nicely fitting 
the Li-mass trend. Models at [Fe/H] - -2.3 including diffusion 
and T6.0 (short dashed line), T6.09 (dotted Une) and T6.25 (solid 
line) turbulence (Richard et al. 2005) are shown. The models 
were rescaled to an initial Au-2.64 (long dashed line) and by 
AM = +0.05 Mq. 
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Table 1. Parameters and Li abundances of our metal-poor stars 



Star 


EW+cr 


ref. 


E(B-V)+o- 


ref. 


Teff 


logg 


[Fe/H] 


mass+cr 


4™ 






(mA) 




(mag) 




(K) 


(dex) 


(dex) 


(Mo) 


(dex) 


(dex) 


BD+01 3597 


24.9±0.5 


A06 


0.045+0.005 


NaD 


6289 


4.04 


-1.87 


0.784+0.018 


2.22 


2.18 


BD+02 3375 


33.3±1.0 


MR04+S07 


0.034+0.010 


NaD 


6163 


4.13 


-2.24 


0.748+0.010 


2.27 


2.21 


BD+02 4651 


31.6±0.8 


UVES archive 


0.038+0.004 


NaD 


6349 


3.79 


-1.78 


0.909+0.062 


2.42 


2.35 


BD+03 0740 


20.8±0.5 


A06 


0.022+0.003 


NaD 


6419 


3.97 


-2.71 


0.781+0.021 


2.21 


2.17 


BD+09 0352 


34.0±1.5 


MR04 


0.012+0.004 


maps 


6145 


4.32 


-2.05 


0.743+0.015 


2.28 


2.21 


BD+09 2190 


16.6±0.5 


A06 


0.015+0.003 


NaD 


6479 


4.01 


-2.64 


0.787+0.030 


2.18 


2.13 


BD+110468 


26.0±1.2 


B05 


0.004+0.004 


maps 


5739 


4.58 


-1.55 


0.693+0.015 


1.85 


1.82 


BD+17 4708 


27.6±0.5 


A06 


0.010+0.002 


NaD 


6127 


4.04 


-1.59 


0.783+0.014 


2.19 


2.14 


BD+24 1676 


24.2±1.1 


UVES archive 


0.013+0.004 


NaD 


6387 


3.84 


-2.54 


0.840+0.071 


2.27 


2.21 


BD+26 4251 


35.5±0.4 


UVES (NSIO) 


0.007+0.003 


NaD 


6131 


4.39 


-1.27 


0.813+0.014 


2.33 


2.26 


BD+29 2091 


37.1±1.0 


UVES archive 


0.004+0.003 


NaD 


5974 


4.58 


-1.99 


0.720+0.014 


2 21 


2.17 


BD+34 2476 


23.8±1.0 


MR04+B05+S07 


0.006+0.005 


maps 


6416 


3.95 


-2.07 


0.805+0.028 


2.30 


2.23 


BD+36 2165 


33.1±0.9 


FIES/NOT 


0.002+0.002 


NaD 


6315 


4.28 


-1.38 


0.846+0.011 


2.42 


2.35 


BD+38 4955 


4.9+0.5 


MR04 


0.000+0.010 


nearby 


5265 


4.60 


-2.59 


0.578+0.010 


0.75 


0.74 


BD+42 2667 


31.7±0.5 


HIRES archive 


0.003+0.001 


NaD 


6148 


4.26 


-1.40 


0.792+0.010 


2 27 


2 21 


BD+42 3607 


39.5±0.5 


HIRES archive 


0.018+0.003 


NaD 


6021 


4.59 


-2.06 


0.733+0.015 


2.29 


2.22 


BD+51 1696 


24.0±0.6 


FIES/NOT 


0.000+0.002 


NaD+nearby 


5722 


4.59 


-1.40 


0.699+0.010 


1.80 


1.77 


BD-04 3208 


23.9±0.5 


A06 


0.009+0.003 


NaD 


6491 


3.98 


-2.29 


0.813+0.021 


2.38 


2.30 


BD- 10 0388 


29.1±0.5 


A06 


0.009+0.003 


maps 


6260 


3.98 


-2.32 


0.785+0.018 


2.27 


2.21 


BD-13 3442 


20.4±0.5 


A06 


0.011+0.002 


NaD 


6435 


3.93 


-2.76 


0.801+0.040 


2.22 


2.18 


CD-24 17504 


18.2±0.5 


HIRES 


0.020+0.005 


NaD 


6451 


4.13 


-3.34 


0.765+0.019 


2.18 


2.14 


CD-30 18140 


27.6±0.5 


A06 


0.012+0.003 


NaD 


6373 


4.13 


-1.90 


0.793+0.010 


2.35 


2.27 


CD-33 01173 


16.0±0.5 


A06 


0.005+0.005 


NaD 


6536 


4.29 


-3.01 


0.786+0.016 


2.18 


2.14 


CD-33 03337 


39.0±0.5 


A06 


0.007+0.006 


NaD 


6001 


4.00 


-1.33 


0.792+0.010 


2.30 


2.24 


CD-35 14849 


28.4±0.5 


A06 


0.002+0.002 


NaD 


6396 


4.22 


-2.35 


0.765+0.010 


2.37 


2.29 


CD-48 02445 


25.6±0.5 


A06 


0.015+0.008 


maps 


6453 


4.25 


-1.93 


0.813+0.016 


2.38 


2.30 


GO 11-044 


33.0±1.2 


B05 


0.006+0.006 


maps 


6304 


4.30 


-2.08 


0.771+0.015 


2.38 


2.30 


G024-003 


28.1±1.3 


UVES archive 


0.014+0.004 


NaD 


6118 


4.27 


-1.59 


0.767+0.012 


2.19 


2.15 


G053-041 


26.3±0.5 


FIES/NOT 


0.009+0.002 


NaD 


6049 


4.31 


-1.24 


0.784+0.012 


2.14 


2.09 


G064-012 


23.7±0.5 


HIRES (AMOS) 


0.003+0.002 


NaD 


6463 


4.17 


-3.26 


0.763+0.010 


2.32 


2.25 


G064-037 


16.4±0.5 


HIRES (AMOS) 


0.012+0.003 


NaD 


6583 


4.20 


-3.17 


0.786+0.012 


2.21 


2.17 


G075-031 


37.8±0.5 


A06 


0.005+0.001 


NaD 


6157 


4.19 


-1.01 


0.883+0.015 


2.40 


2.33 


Gl 14-042 


25.2±1.2 


UVES (NSIO) 


0.012+0.002 


NaD 


5848 


4.40 


-1.09 


0.762+0.015 


1.91 


1.89 


G 192-043 


29.6+1.4 


FIES/NOT 


0.012+0.003 


NaD 


6298 


4.39 


-1.39 


0.844+0.017 


2.35 


2.27 


HD003567 


38.3±0.5 


A06 


0.004+0.003 


NaD 


6177 


4.14 


-1.17 


0.846+0.010 


2.41 


2.34 


HDO 16031 


28.4±0.6 


UVES archive 


0.003+0.003 


NaD 


6286 


4.17 


-1.74 


0.789+0.011 


2.29 


2.23 


HDO 19445 


34.9±0.5 


A06 + HIRES 


0.000+0.000 


NaD+nearby 


6136 


4.45 


-2.02 


0.739+0.011 


2.29 


2.22 


HD024289 


47.6+1.7 


MR04 


0.015+0.006 


maps 


5927 


3.73 


-2.08 


0.798+0.054 


2.31 


2.24 


HD029907 


26.0+1.3 


MR04 


0.000+0.000 


nearby 


5531 


4.63 


-1.58 


0.652+0.010 


1.69 


1.67 


HD031128 


30.9+1.1 


UVES archive 


0.000+0.000 


nearby 


6092 


4.51 


-1.52 


0.771+0.011 


2.22 


2.17 


HD034328 


35.8+1.0 


UVES archive 


0.000+0.000 


NaD+nearby 


6056 


4.50 


-1.66 


0.757+0.012 


2.27 


2.20 


HD059392 


39.2+0.5 


A06 


0.006+0.003 


NaD 


6045 


3.87 


-1.62 


0.792+0.012 


2.32 


2.25 


HD064090 


12.1+1.0 


MR04 


0.000+0.000 


nearby 


5465 


4.61 


-1.68 


0.642+0.010 


1.27 


1.26 


HD074000 


24.0+1.0 


UVES arc+MR04 


0.003+0.002 


NaD 


6362 


4.12 


-2.01 


0.779+0.010 


2.25 


2.20 


HD084937 


25.3+0.5 


HIRES (AMOS) 


0.005+0.002 


NaD 


6408 


3.93 


-2.11 


0.786+0.010 


2.32 


2.26 


HD094028 


36.9+0.9 


UVES archive 


0.000+0.000 


NaD+nearby 


6111 


4.36 


-1.45 


0.777+0.010 


2.33 


2.26 


HD 102200 


33.0+0.5 


A06 


0.005+0.004 


NaD 


6155 


4.20 


-1.24 


0.822+0.010 


2.30 


2.24 


HD108177 


31.2+0.9 


UVES archive 


0.003+0.002 


NaD 


6334 


4.41 


-1.65 


0.812+0.010 


2.40 


2.31 


HD 11 6064 


29.7+0.6 


UVES archive 


0.000+0.003 


NaD+nearby 


5896 


4.41 


-1.89 


0.702+0.010 


2.03 


1.97 


HD122196 


40.7+0.6 


UVES archive 


0.004+0.002 


NaD 


5986 


3.73 


-1.81 


0.782+0.010 


2.29 


2.23 


HD126681 


14.0+0.4 


UVES (NSIO) 


0.000+0.001 


NaD+nearby 


5639 


4.57 


-1.17 


0.714+0.013 


1.48 


1.48 


HD 132475 


56.6+0.5 


UVES archive 


0.007+0.003 


NaD 


5746 


3.78 


-1.52 


0.808+0.015 


2.25 


2.20 


HD 140283 


47.5+0.5 


A06+AM08 


0.000+0.002 


NaD+nearby 


5777 


3.62 


-2.40 


0.774+0.010 


2.17 


2.13 


HD160617 


40.5+0.5 


A06 


0.005+0.004 


NaD 


6048 


3.73 


-1.78 


0.820+0.023 


2.33 


2.27 


HD163810 


21.6+0.7 


UVES (NSIO) 


0.024+0.006 


NaD 


5602 


4.59 


-1.29 


0.692+0.013 


1.64 


1.63 


HD166913 


38.7+0.6 


UVES archive 


0.000+0.000 


NaD+nearby 


6268 


4.29 


-1.57 


0.799+0.010 


2.46 


2.39 


HD181743 


38.0+1.4 


MR04 


0.007+0.007 


M06 


6151 


4.48 


-1.78 


0.760+0.013 


2.36 


2.28 


HD189558 


58.2+0.7 


UVES archive 


0.000+0.004 


NaD+nearby 


5765 


3.78 


-1.12 


0.863+0.014 


2.29 


2.24 


HD 193901 


28.0+0.6 


UVES archive 


0.001+0.001 


NaD 


5915 


4.51 


-1.07 


0.785+0.014 


2.05 


2.00 


HD 194598 


31.5+0.5 


HIRES 


0.000+0.000 


NaD+nearby 


6118 


4.37 


-1.13 


0.828+0.010 


2.26 


2.20 


HD 199289 


19.7+0.8 


UVES archive 


0.000+0.000 


NaD+nearby 


5975 


4.35 


-1.02 


0.798+0.010 


1.93 


1.90 


HD201891 


24.3+0.5 


HIRES 


0.000+0.000 


nearby 


5947 


4.31 


-1.05 


0.788+0.010 


2.01 


1.96 


HD205650 


15.6+0.5 


UVES (NSIO) 


0.000+0.000 


NaD+nearby 


5952 


4.47 


-1.16 


0.779+0.013 


1.82 


1.78 


HD2 13657 


30.0+0.5 


A06 


0.002+0.002 


NaD 


6299 


3.90 


-1.96 


0.795+0.023 


2.33 


2.26 


HD2 18502 


28.6+0.8 


UVES archive 


0.000+0.005 


NaD+nearby 


6298 


3.96 


-1.76 


0.791+0.010 


2.31 


2.24 


HD219617 


37.8+0.8 


UVES arc+MR04 


0.001+0.002 


NaD 


6146 


4.10 


-1.48 


0.790+0.012 


2.37 


2.30 


HD233511 


31.9+1.5 


FIES/NOT 


0.002+0.002 


NaD 


6127 


4.29 


-1.63 


0.759+0.010 


2.25 


2.19 


HD241253 


19.8+0.4 


UVES (NSIO) 


0.001+0.001 


NaD 


6023 


4.37 


-1.05 


0.818+0.010 


1.98 


1.93 


HD250792 


16.0+0.7 


FIES/NOT 


0.002+0.002 


NaD+nearby 


5568 


4.40 


-1.01 


0.715+0.010 


1.50 


1.50 


HD284248 


32.0+1.0 


UVES archive 


0.020+0.005 


NaD 


6285 


4.29 


-1.56 


0.810+0.010 


2.37 


2.29 


HD298986 


30.4+0.5 


A06 


0.004+0.002 


NaD 


6223 


4.26 


-1.34 


0.825+0.010 


2.30 


2.23 


HD338529 


23.5+0.5 


A06+HIRES 


0.008+0.002 


NaD 


6426 


3.89 


-2.29 


0.795+0.018 


2.30 


2.23 


LP0635-0014 


22.3+0.5 


HIRES 


0.060+0.010 


NaD 


6516 


3.98 


-2.50 


0.808+0.041 


2.35 


2.28 


LP08 15-0043 


18.9+0.5 


A06 


0.024+0.015 


M06 


6534 


4.17 


-2.78 


0.781+0.010 


2.26 


2.20 



